Generalized model-independent approach to intrinsic decoherence 
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A formalism is presented to express decoherence both in the markovian and nonmarkovian regimes 
and both dissipative and nondissipative in isolated systems. The main physical hypothesis, already 
contained in the literature [Q, ^, amounts to allowing some internal parameters of the system to 
evolve in a random fashion. This formalism may also be applicable to open quantum systems. 
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Decoherence is at the heart of the destruction of quan- 
tum superpositions in most physical systems. Its role is 
basic in questions such as the quantum-classical transi- 
tion the realization of quantum computation or the 
quantum processing of information Q. From a quan- 
titative point of view the effect of decoherence is com- 
monly described by the action of the environment upon 
the physical system of interest |^. This view in- 
evitably depends on the chosen modelization of this envi- 
ronment. Recently |^ (see also 0) a model-independent 
approach to decoherence has been initiated based on the 
hypothesis of randomness of some of the internal param- 
eters of the system, thus reducing the role of the en- 
vironment. Here we present a triple generalization of 
this model-independent approach which allows us to find 
more general master equations both in the markovian 
and in the nonmarkovian regimes] 17 and both dissipa- 
tive and nondissipative. 

Randomness in an isolated system is introduced by 
claiming that the evolution time is a stochastic parame- 
ter I?! (previous references in this respect are fl], ||). 
This physical hypothesis is implemented by noticing that 
the experimental observations are described by the time- 
averaged density operator, as expressed by the relation 
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dt'P{t,t')p{t') 
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where P{t, t') denotes the probability density associated 
with the random time t' and p{t') denotes the usual quan- 
tum density operator. To find an expression for P{t,t') 
two conditions are imposed |^ : i) p{t) must be a density 
operator itself; ii) the evolution of p{t) must be marko- 
vian. Under these assumptions and using the integral 
definition of the Gamma function, one arrives at a fixed 
probability density P{t, t') (up to two parameters) and 
thus at a fixed time-averaged density operator p{t) . Here 
we drop condition ii), thus comprising nonmarkovian evo- 
lutions and moreover using different tools we also allow 



for distinct markovian situations. The nonmarkovian 
regime is interesting in its own since it is the rigorous 
type of dynamics (prior to approximations) followed by 
a physical system coupled to an environment |l^ . Fur- 
thermore the experimental advances in mesoscopic scales 
(for instance, in creating atom lasers with Bose-Einstein 
condensates [^) also demands a better understanding 
of this regime. The generality obtained also allows us to 
embrace dissipative decoherence in contrast to 0, 0] . 

The idea is to express the evolution operator Ujt) us- 
ing the spectral decomposition theorem for unitary oper- 
ators [T^ and then introduce the randomness in the time 
evolution. Let us recall then that 
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where Eg is the spectral family on [~tt, tt] associated to 
U{t). The randomness hypothesis can be introduced in 
different ways. The most nai've way is just to substitute t 
in expression (|^) for a random variable Xt which expresses 
the stochastic nature of the evolution. The properties 
of Xt (distribution, moments, ... ) will depend on the 
characteristics of the physical system under study. To 
grasp this generality it is convenient to consider xt as the 
solution to an Ito stochastic differential equation |]l3| : 



dxt = Xt)dt + a{t, Xt)dBt 



(3) 



where 6(-, ■) and a{-, ■) are real functions of two real vari- 
ables and Bt denotes standard real brownian motion. 
These two functions encode the information of the spe- 
cific physical system. As a first example we may consider 
K^jXt) = 1 and a{t,xt) = cr{t) (with clear abuse of no- 
tation). This allows us to write Xt = t + Jq a{s)dBs, 
hence the expectation value of Xt is the standard de- 
terministic time t and its deviation from it is given by 
A(t) = j^a^{s)ds. This can be considered a weak- 
coupling- like limit. With these hypotheses, we may con- 
struct a random evolution operator given by 
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J — TT 
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Following [g| the density operator describing experi- 
mental results will be given by 

m = E[p{t)] - E[C/,and(Op(0)t/Ld(i)] (5) 

where E denotes expectation value with respect to the 
probability distribution of xt- Using Ito calculus [p^ , 
the expectation value E[e-'^'^'^'^"> "^^b,] = e-^^^'"')' 
can be readily found and hence one immediately arrives 
at the following expression for p{t) 

=e-**^«-^^«[p(0)] (6) 

where Ch = [H, ■] is the commutator with the hamilto- 
nian H. Note that this equation automatically satisfies 
condition i) above, i.e. p{t) is a density operator (it is 
selfadjoint, it has trace one and it is positive) provided 
the initial p(0) is a density operator too. This is a first 
generalization of the results obtained in If condi- 
tion ii) is imposed, then the only choice for ait) will be 
a{t) — 7^/^, i.e. a constant. The master equation is the 
well-known phase-destroying master equation 

m = -^[H,m]-l[H,[H,m]] m 

Note that no approximation has been assumed to ar- 
rive at this equation, in contrast to [Q; ^. The constant 
7 is thus under this formalism a measure of the devia- 
tion from strictly determinist time evolution. Apart from 
the evident freedom in the choice of adequately physically 
motivated choices of functions &(•, •) and cr(-, •) (which e.g. 
may drive us to time-dependent 7(t)), a higher general- 
ization can be achieved if the previous scheme is slightly 
changed by substituting the whole phase of the exponen- 
tial in (|) for a 0-dcpcndcnt random variable xt (^) ■ This 
allows us to write the random evolution operator instead 
as 



C/,a„d(t) - / e-'^^^'UEe (8) 



The expectation value is found using standard tech- 
niques. Define the stochastic process Xt{0,9') by the 
SDE dXt{e,e') = a{t-e)dBt{9) - a{t;9')dBt{9'); express 
exp(— iXj(6', 6*')) as a power series; apply the Ito formula 
(cf. H) to f{Xt) ifiz) = z"); define f3„{t;9,9') = 
E[Xl'(9,9')] to find 



The former scheme is a particular case of this latter 
when Xt(^) — ^Xt- The physical meaning behind this 
formalism is clear: the random time evolution does not 
affect on an equal footing to all energy levels of the sys- 
tem. Following then the same idea as above, we consider 
Xt{t) as a solution to an Ito SDE, now 6'-dependent 



dxtie) ^ b{t,xt{0);9)dt + a{t,xt{ey,9)dBt{9) (9) 



i.e., now we have a SDE for each energy level. This 
must be supplemented with correlation properties among 
the brownian motions at different 9^s. These will be ex- 
pressed through the relation 



dBt{9)dBt{9')=g{t-9,9')dt (10) 



where g{t] 9, 9') is basically for each t a covariance func- 
tion expressing the correlations among different Bt{9ys. 
From a physical standpoint it is clear that g{t; 9, 9') = 
g{t] 9' , 9) and since we are dealing with standard brown- 
ian motions \g{t\9,9')\ < 1 for all t,9,9'. 

It is difficult to proceed forward without making some 
assumptions. One may on one hand discuss on physical 
grounds particular choices for the different functions in- 
volved and find the solution or on the other hand look 
for general analytical properties which drives us to a suf- 
ficiently global situation. As before, as a first example, 
let us assign simple choices to h{t,xt{0)]9) = h{9) and 
a{t,xt{0)',9) — a{t\9). The former corresponds to the 
fact that the unitary part of the dynamics is driven by a 
renormalized hamiltonian [Q, as the following calcula- 
tion shows; whereas the second states that the departure 
from deterministic evolution is caused by an uncorrelated 
noise with different intensity at different times for each 
9. Let us include some details. We have to evaluate the 
expression: 



(11) 

I 

f32n{t;9,9') = ^—>-r{t;9,9') (12a) 
I32n+iit;9,9') = (12b) 

where X{t;9,9') = j^ri{s;9,9')ds and we have de- 



-p{t)= f f e-'(Mfl)-Me'))*E[e-</o-(«;e)d8,W-/X.;e')<iB.(e'))]^^^^(0)rf^^ 
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fined r){t;9,9') = a^{t;9)g{t;0,e) + a^{t;e')g{t;e' ,9') ^ 
2a{t; 9)a{t; 9')g{t; 9, 9'). Finally resum the series and 
then ( pi] ) reduces to 

J —T: J —TT 

(13) 

We cannot proceed further without specifying a{t] 9) 
and g{t] 9, 9'). To compare with the previous scheme let 
us define (T(t; 9) = 7^/^6' and assume a correlation func- 
tion git;9,9') = exp{-T'^{9 - O'f). Here r'^ roughly 
measures the distance of correlation in some arbitrary 
units. Then X{t;9,9') = -/t{9^ + 9'^ - 299' e"'^^'>''>'^^) 
and the master equation reads 

p{t) - -i[h{H)rp{t)]- 

~ l{H'-p{t) + -p{t)H^ - 2i/e--'^«[p(i)]i/)(14) 

This is a CP markovian master equation [p^ different 
from (^, though it reduces to it in the limit r — > 0, 
i.e. in the case of infinite-range correlations (all en- 
ergy levels equally affected by the random evolution). 
In this way the formalism depicted here generalizes the 
approach in even when restricting ourselves to the 
markovian regime. Indeed, whenever we may write 
A(t; 9, 9') = tX{9, 9') (which amounts to having (j{t; 9, 9') 
and g{t;9,9') time-independent) we arrive at a marko- 
vian evolution, since 

m = r r e*[-^(''(^)-"(«'))-^(«'«')idij,p(o)dii;,, 

J —TT J —IT 

= e-"^M^f)+*i'«[p(0)] (15) 

where Ch{H) and Vh denote superoperators. 

The previous example depicts a typical situation of 



nondissipative decoherence. In this language and re- 
stricting to the previous case (eq. (^), this is due to 
the fact that X{t,9,9) = 0, as an immediate calculation 
of Tr{p(t)H) shows. Working in the energy eigenstates 
basis, (H) leads to p{t;E,E') = exp{-it{hiE)-h{E'))- 
X{t■, E, E')/2)p{Q- E, E'), hence 



Tiip{t)H) = I Ee-''^'''^'^%{0;E,E)dp{E) = 

Ja{H) 

= Tr{p{0)H) (16) 

where (t{H) denotes the spectrum of H and fi{E) denotes 
the measure over this spectrum (a sum if it is discrete and 
an integral if it is continuous) and where A(t, 9,9) = has 
been assumed. From (|l^) it is clear that if this condition 
is not satisfied, we have a dissipative decoherence pro- 
cess. This is a third generalization of previous works. 
Some comments must be made about this formalism. 
Though it has been presented as an adequate tool to 
express decoherence in isolated systems, from a mathe- 
matical standpoint it can be viewed as an alternative de- 
scription of environment-induced decoherence provided 
the stochastic parameter xt is properly connected to the 
surrounding reservoir properties and the environment- 
system interaction characteristics. But in our opinion 
the main benefit from this language is to bring a bit 
closer the decoherence process and the different models of 
stochastic state reduction (see for a review of their 
common backbone). However this connection demands 
further mathematical work which is in progress. 
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